Multiscale behaviour of volatility autocorrelations in a financial market 
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We perform a scaling analysis on NYSE daily returns. We show that volatility correlations are 
power-laws on a time range from one day to one year and, more important, that they exhibit a 
multiscale behaviour. 
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It is well known that stock market returns are uncor- 
related on lags larger than a single day, in agreement 
with the hypothesis of efficient market. On the contrary, 
absolute returns have memory for longer times; this phe- 
nomenon is known in financial literature as clustering 
of volatility. In ARCH-GARCH models volatility 
memory is longer than a single time step but it decays ex- 
ponentially while empirical evidence is for hyperbolic cor- 
relations (4^. In this paper, we perform a scaling analy- 
sis of the standard deviation of a new class of observables, 
the generalized cumulative absolute returns. This analy- 
sis clearly shows that volatility correlations are power- 
laws on a time range from one day to one year and, more 
important, that the exponent is not unique. This kind of 
multiscale behaviour is known to be relevant in the theory 
of dynamical systems, of fully developed turbulence and 
in the statistical mechanics of disordered systems (see 
fiof for a review) while it is a new concept for financial 
modeling. 

We consider the daily New York Stock Exchange 
(NYSE) index, from January 1966 to June 1998, for a 
total of N = 8180 working days. The quantity we con- 
sider is the (de-meaned) daily return, defined as 
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where St is the index value at time t ranging from 1 to 
N, and (•} is the average over the whole sequence. The 
underlying daily volatility at is not directly observable, 
but it is indirectly defined by r t = er t ?7 t . It is assumed 
that the r\t are identically distributed random variables 
with vanishing average and unitary variance. The usual 
choice for the distribution of the r\t is the normal Gaus- 
sian. The observables directly related to the volatilities 
are the absolute returns \r t \. 

As pointed out by several authors 0-|l3|, the distri- 
bution of returns is leptokurtic. In p2| , it was firstly 
proposed a symmetric Levy stable distribution and more 
recently in [ [l3[ it has been provided strong evidence for 
this fact. More precisely, in Jl3] it is shown that the dis- 
tribution is Levy stable for high frequency returns except 
for tails, which are approximately exponential. The esti- 
mation is that the shape of a Gaussian is recovered only 
on longer scales, typically one month. 

Our analysis is on low frequency data, and first of all 
we want to verify that anomalous Levy scaling is not 



effective in a range of time from one day to one year. 
We consider the cumulative returns <f>t(L), defined as the 
sum of L successive returns ft, ... , rt+L-i, divided by L. 
Using NYSE data one can define N/L not overlapping 
variables of this type and compute the standard deviation 
o~(L). The standard deviation is independently computed 
for L ranging from 1 to 250 (one year). Larger value of L 
would imply insufficient statistics. Assuming that r t are 
uncorrected (or short range correlated), it follows that 
o~(L) has a power-law behaviour with exponent 0.5 for 
large L, i.e. cr(L) ~ with (3 — 0.5. The exponent for 
the NYSE index turns out to be about 0.49 (see fig. 1 and 
also see 14 ) , according to the hypothesis of uncorrelated 
returns. This value of the exponent also ensures that 
Levy scaling is not effective in this range of time. 
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FIG. 1. Standard deviation S(L,y) of the generalized cu- 
mulative absolute returns (2) as a function of L on log-log 
scales for 7 = 1 (crosses) and 7 = 1.5 (slanting crosses), 
compared with the standard deviation of the cumulative re- 
turns (circles). The exponents of the best fit straight lines 
(dashed lines) are, respectively, 13(1) ~ 0.19, /3(1.5) ~ 0.27 
and ~ 0.49. 

Let us introduce the generalized cumulative absolute 
returns defined as 
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where 7 is a real exponent and again, these quantities are 
not overlapping. If the |r t | 7 are uncorrelated, one should 
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find that the standard deviation 6(L, 7) has a power-law 
behaviour with exponent 0.5. 

On the contrary, a power-law autocorrelation function 
with exponent a( 7 ) < 1 (N> t+L | 7 )-<|r t | 7 ><|r t+L | 7 ) ~ 
L~ a ^\ would imply that 5(L, 7) is a power-law with ex- 
ponent /3( 7 ) = a(-f)/2. For autocorrelations with expo- 
nent a (7) > 1 we would no detect anomalous scaling for 
the standard deviation (/3( 7 ) = 0.5). 

Our numerical analysis on the NYSE index shows very 
sharply that S(L, 7) has an anomalous power-law be- 
haviour in the range from one day to one year (L — 250). 
For example, for 7 = 1 we find (3(1) ~ 0.19, while for 
7 = 1.5, /3(1.5) ~ 0.27 (sec fig. 1). For larger L the 
statistics becomes insufficient. 

The crucial result is that /3( 7 ) is a not constant func- 
tion of 7 in the range —0.5 < 7 < +4 (see fig. 2), show- 
ing the presence of different scales. The interpretation is 
that different 7 select different typical fluctuation sizes, 
any of them being power-law correlated with a different 
exponent. 
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autocorrelations would not have provided an analogous 
clear evidence for multiscale power-law behaviour, since 
the data show a wide spread compatible with different 
scaling hypothesis. 
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FIG. 3. Autocorrelation function of \rt\ (crosses) as a func- 
tion of the correlation length L, compared with the autocor- 
relation function of r t (circles). The data are in agreement 
with a power-law with exponent 2/3(1) ~ 0.38 in the first case, 
and absence of correlations in the second. In both cases the 
scale is fixed by autocorrelations equal to 1 at L — 1. 

In conclusion, we have found that scaling of variance 
of the generalized cumulative absolute returns implies 
multiscaling power-law correlations in financial indices. 
This result clearly suggests that models with exponential 
correlations, like ARCH-GARCH, are inadequate to de- 
scribe the dynamics of financial markets, and they should 
be implemented to account for the coexistence of long 
memory with different scales. 



FIG. 2. Scaling exponent f3(-y) of the standard deviation 
5(L, 7), where the bars represent the errors over the best 
fits. An anomalous scaling (/3 < 0.5) is shown in the range 
-0.5 < 7 < +4. 

The longest correlation is for 7 = 0.15 (/?(0.15) ~ 
0.15). The case 7 = corresponds to cumulative log- 
arithm of absolute returns. 

In the region 7 > 4 the averages are dominated by 
only few events, corresponding to very large returns, and, 
therefore, the statistics becomes insufficient. 

The anomalous power-law scaling can be directly 
tested against the plot of autocorrelations. For instance, 
the autocorrelations of r t and of \r t \ are plotted in fig. 3 
as a function of the correlation length L. Notice that 
the full line, which is in a good agreement with the 
data, is not a plot but it is a power-law whose exponent 
2/3(1) ~ 0.38 is obtained by the previous scaling analysis 
of the variance. The autocorrelations for the return, as 
expected, vanish except for the first step (L = 1). 

It should be also noticed that a direct analysis of the 
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